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Extensive numerical results for electromagnetic scattering; from spheres are presented 
in graphical form. To simulate various types of coatings, the surface impedance of the 
sphere is a specified complex function. It is shown conclusively that the forward scattering 
lobe is almost independent of the characteristics of the coating. On the other hand, the 
backscattering and side scattering are profoundly modified by changes in the effective 
impedance of the coating. 

1. Introduction 

Scattering of electromagnetic waves from spheres is a subject which has permeated the 
literature for over half a century. An exhaustive survey of the subject is found in a highly 
recommended book by Van de Hulst [1957.] An important collection of informative papers 
reporting recent advances is also available [Kerker, 1964]. 

Despite the numerous papers treating both theoretical and experimental aspects, a number 
of unresolved questions remain. For example, the effect of coating conducting spheres with 
absorbent material is not entirely understood. Also, the related question of how an inductive- 
or capacitive-type surface layer modifies the scattering properties has not been answered 
completely. It is the purpose of this paper to present calculations for scattering from spheres 
under various surface impedance boundary conditions. This relatively simple model has the 
virtue of displaying many important features in an economical manner. This approach has 
been used previously by various investigators [Wait, 1956; Hiatt et al., 1960a, 1960b; 
Garbacz, 1964] who have treated closely related problems. 

2. Formulation 

The spherical model used as a basis for calculation is illustrated in figure 1 with reference 
to a Cartesian coordinate system (x, y, z). A plane wave is incident from the negative z 
direction and polarized with the electric field E Q in the x direction. With reference to a spherical 
coordinate system (/•, 6, (/>) also centered at the sphere, the scattered field components may be 
written 

E- e ==^e-*> (cos 4>)P(e), (1) 



and 
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E%=^e-«" (sin 4>)Q(d), (2) 




INCIDENT PLANE WAVE 



Figure 1. Geometry for the spherical scatterer*. 



where k= 2 ir /wavelength and where P and Q are functions of only. The latter are defined by 



«»-a^j[*i«<-«+«-. a £J fi } 



and 



<2W=S 



2w+ 



Z=in(n+ 



i r R Pi ( 



( cos9 ) + C > ^P' (cose) 



] 



(3) 
(4)^ 



where the coefficients B n and C w are a function of ka, the circumference of the sphere in wave- 
lengths, and the surface impedance Z of the sphere. The explicit definition of these coefficients 
B n and C n is given in the appendix where the derivation is outlined. The function P\ (cos 0) is 
the associated Legendre function of order n and index 1 . 

As indicated above, the solution for the present problem is expressed in terms of the surface 
impedance Z which requires that the total tangential fields at r=a are related as follows: 



Ee=-ZH<, and E*=ZH e . 



(5) 



These are often termed the Leontovich boundary conditions when applied to the interface 
between free space and a highly conducting medium [e.g., ch. 2 in Wait, 1962]. Actually, they 
are a valid representation for other type surfaces such as a dielectric coated conductor or a 
finely corrugated surface, in which cases Z is highly inductive. In the general case, such as 
for a radially inhomogeneous sphere, it is necessary to recognize that surface impedances are 
dependent on mode number. 

For the present discussion, it is assumed that the boundary may be characterized by a< 
constant surface impedance Z. This is conveniently normalized by writing 



Z/ 7) = iG, 
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(6) 



where 7j =l20ir is the intrinsic impedance of free space and 6r is a dimensionless factor. For a 
purely inductive surface, G is positive real, while for a capacitive surface, G is negative real. 
If the sphere is composed of a homogeneous, highly conducting material, 

0"^|0|> (7) 

where \G\ 2 = e ufa g in terms of the conductivity a g of the sphere, the dielectric constant of free 
space e , and the angular frequency co. The restrictions inherent in the use of the surface 
impedance require that |#| 2 <<1 [Wait, 1962] in (7). 

In the general case of a surface impedance boundary, it is desirable to write 



-or 



Z=|Z|e* (8) 

G=-i\G\e l \ (9) 



where $ is a phase angle. For the inductive boundary $=90°, for the capacitive boundary 
<t>=— 90°, whereas, for the homogeneous conductor $=45°. It might be mentioned that witli 
various forms of concentrically stratified conducting spheres <£ may be anywhere in the range 
to 90°. Furthermore, certain types of corrugated surfaces loaded with lossy material could 
give rise to effective surface impedances which have phase angles also in this range [Wait, 1962]. 
To present numerical results in an economical manner, two bistatic radar cross sections 
<7 t (6) and <r h (9) are introduced which are for Emplane and //-plane scattering. Explicitly, these 
are defined by: 



and 



<*h\P)= W~2 • (11) 

■&Q J = 90° 



Using (1) and (2), and dividing by the geometrical cross section, it is seen that (10) and (11) 
may be written 

°3=w^ m > (12) 

.and 

if-TK> w - <13) 

A computer program was devised for the functions P(6) and Q(6). The spherical Bessel 
functions which are involved in the coefficients B n and C n were computed by a method due 
essentially to Miller [1950]. To secure convergence of the series expansions for P(6) and Q(6), 
at least 2 ka terms were used in the summation. For the graphical accuracy shown in the 
figures, this was more stringent than necessary. 

3. Description of Results 

The numerical results for the two bistatic cross sections are plotted as a function of 6 for 
a wide variety of situations. The basic parameters are the circumference, ka, of the sphere 
in wavelengths and the complex value of G. 
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Figure 2a. 




Figure 2b. 



In figures 2a and 2b, the cross sections are plotted for ka values ranging from 3 to 20 under 
the condition 67=0 (i.e., perfectly conducting sphere). These particular curves may be com- 
pared directly with similar plots prepared by Nelson Logan 1 and colleagues. Where the ka 
values are common the data are identical within the graphical accuracy. This would appear 
to be a good check on the computer program for this special case. 

It is evident from the curves in figures 2a and 2b that the forward scattering (i.e., cor- 
responding to 0^0) is greatly enhanced over the backscattering. In fact, the (normalized) 
forward scattering cross sections both become approximately equal to (ka) 2 while the (normal- 
ized) backscatter is approximately unit. Thus, in the forward direction, the sphere is behaving 
as a "black disk," particularly for the larger ka values. On the basis of HuygeiVs principle, 
it follows from Van de Hulst [1957] that for small values of 6, 



<r e (0)^<r h (0)c±(ka) 2 IT e -M*°°**+v*toV*inedxdy, 



(14} 



where the integration is over the area of the "black disk." In the case of the sphere or its- 
equivalent circular black disk, the integrations are readily carried out to give 



e (e)a h (d) 



7r(r 



7r<r L (ka sm B) J 



(15) 



For pure forward scattering (i.e., 0~ K)), the square bracketed term involving the J x Bessel 
function approaches unity which is a broad maximum. It decreases smoothly to zero at 6=6 , 
when Ji(ka sin d )=Q or when ka sin O =3.832. For ka = 10, 15, and 20, the corresponding 
values of 6 are 22.5°, 14.8°, and 11.0°. These correspond quite well to the first minima in the 
calculated curves (fig. 2b) for a e (6) and a h (6) in the forward direction. The agreement is better 
for the larger ka values, as expected. 



i These curves appear in the appendix to a monograph by King and Wu [1959]. 
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Because of the crude physical assumptions underlying the application of Huygen's 
principle, it is not expected that (15) should be used for angles except near forward scattering. 
Furthermore, being a scalar theory, it does not distinguish between the two wave polarizations. 
However, this simple "black disk" model is particularly appropriate for describing the shape 
and width of the Forward main lobe in the scattering pattern. A further demonstration is to 
employ (15) to calculate the values of 0=0 X where <j c {6) and v h {d) drop to 1/2 of the central 
maximum at 6=0. This condition is ka sin ^ = 1.61. For ka=10, 15, and 20, the correspond- 
ing values of 6 { are found to be 9.3°, 6.1°, and 4.6° which, within two significant figures are 
identical to the curves in figure 2b computed from the exact series solutions. 

It is evident that the curves in both figures 2a and 2b are highly oscillatory. This results 
from the well-known effect of the "creeping waves." In other words, the two waves diffracted 
around the sphere interact with one another to form an interference pattern. This effect 
has been discussed for perfectly conducting spheres by Van de Hulst [1957], Belkina and 
Weinstein [1957], Logan and Yee [1961], and many others. It is further remarked here that 
the interference pattern is much more pronounced for .E-plane scattering than it is for S-plane 
scattering. For perfectly conducting spheres, this is also a very well-known effect and is 
primarily a consequence of relatively low attenuation of the creeping wave in the 2£-plane case. 

For perfectly conducting spheres, it is finally remarked that backscattering cross sections, 
for the larger ka values, become very close to the geometrical cross section ira 1 . This is 
particularly evident for the curves in figure 2b, as remarked above. 

The effect of choosing a succession of real positive G values is illustrated in figures 3a to 
5b for the same range of ka values used for the perfectly conducting spheres (i.e., G=0). As 
mentioned before, a purely inductive boundary is described by real positive G. It is very 
evident from these curves that, even for the largest G value chosen, the magnitude and shape 
of the forward -scatter lobe are not changed markedly. Furthermore, the magnitude of the 
backscattering cross section is only slightly modified from the corresponding curves for G=0. 
The only significant changes in the curves in figures 3a to 5b from those in figures 2a and 2b 
are the enhancement of the ripples in the 2?-plane scattering pattern. This is consistent with 
the expected reduction of the attenuation of the creeping wave which results from an 
inductive-type coating [Wait, 1964]. 




Figure 3a. 




Figure 3b. 
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Figure 4a. 
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Figure 4b. 
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Figure 5a. 
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Figure 5b. 



4. Digression on Fresnel Reflection Coefficients 

In discussing the numerical results for backscattering from spheres with impedance 
boundary layers, it is useful to consider first what one may expect on the basis of geometrical 
optics. The concept is extremely simple; it is assumed that the cross sections a e and cr h are 
individually proportional to |i? e | 2 and |/^| 2 where the latter are the square of the Fresnel 
reflection coefficients for S-plane and ZJ-plane scattering, respectively. Thus, within this 
geometrical-optical approximation, 



£s |fl e | 2 and ^ \R h \\ 

TTLl Tib 



(16) 
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The reflection coefficients are defined by [Wait, 1962], 

Re= C+A &ndBh== CA+i' (1?) 

where A = Z/i) n =iGis the aormalized surface impedance and C=cos /3=sin (6/2) where /S is the 
local angle of incidence for the specular ray. 

In terms of the magnitude \G\ and the phase angle $, it is seen that, according to geomet rical 
optics, 

*e_ [sin (g/2) — | g| cos ^] 2 + [ | g| sin *] 2 . 

*a*~[sin (0/2) + \G\ cos <J>] 2 +[|<?| sin *] 2 ' { } 

and 

a^ [\G\ cos * sin (6I2)-If+[\G\ sin $ sin (fl/2)] 2 
*a a ~[\Q\ cos* sin (0/2) + l] 2 +[|£| sin* sin (0/2)] 2 ' 



(19) 



It is rather interesting to note that when G is real (i.e., $ = 90°), the right-hand sides of the 
preceding two equations are unity. This is not really surprising as reflection takes place 
without loss from a purely reactive surface. This is also true, of course, when |6r|=0, which 
corresponds to the perfectly conducting surface. 

5. Further Description of Results and Comparisons With Geometric Optics 

It is now worth while to return and examine the various curves in figures 2a to 5b. It is 
evident that for the larger ka values (i.e., 10, 15, and 20) and for 6 in the range 90° to 180°, 
(T c /iia 2 and cn/tra 2 are within the range 0.6 to 1.5. As expected, of course, geometric optics 
ignores the diffracted or creeping waves. 

To illustrate a much more favorable comparison between the exact theory and geometrical 
optics, we choose a surface impedance which is purely real. Thus, Z=\Z\ and G=—i\G\ so 
that the phase angle <f>=0. For this situation, the scattering patterns are shown in figures 
6a and 6b for the E'-plane and //-plane, respectively. The values of G are indicated on the 
curves. The close agreement between the geometrical-optical calculations and the exact 
theory is clearly evident. It is particularly interesting to note that when |6r| = 1.00, the 
backseat tering according to geometric optics, is identically zero. In other words, the surface 
impedance matches the intrinsic impedance of free space. The exact theory for the sphere 
(with ka=\5) shows that the baekscattering cross section is less .than 10~ 5 . It may also be 
observed in figures 6a and 6b that deep nulls occur whenever |(?|=sin 6/2 in the case of 2?-plane 
scatter and whenever l/|6! r |==siii 6/2 in the case of //-plane scatter. Under these conditions, 
the Fresnel reflections have a zero and the corresponding angle (6/2) is the conventional 
Brewster angle. 

The curves in figures 7a to 9b illustrate the behavior of the cross sections for values of G 
which are near —i. It is seen in figures 7a and 7b that even for a ka value as small as 3, almost 
perfect absorption takes place in the backscatter direction if G=—i. However, very small 
changes in the magnitude of G result in an appreciable backscatter signal with a subsequent 
shift of a deep null to another angle. Similar effects are seen in figures 8a, 8b, 9a, and 9b 
for larger ka values. In all cases, the forward scattering lobe is virtually unchanged even 
when the sphere has no backscatter. This is in agreement with the results of Garbacz [1964] 
for a similar situation. 

In figures 10a to lOe, cross-section curves are shown to illustrate the effect of varying the 
magnitude of G over fairly wide limits. Both E- and //-plane cases are shown for ka values 
of five and fifteen. In figure 10a, the value of G is +0.50, which corresponds to a highly 
inductive surface. Here, there is evidence of strong standing waves for the 22-plane pattern 
for the whole range of 6. As explained before, this results from the low attenuation of the 
creeping waves for the TM (transverse magnetic) waves. On the other hand, the standing- 
waves for the //-plane pattern are only slightly enhanced by the inductive coating except for 
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Figure 6b. 



directions near backscatter. As seen from figure 10b, the effects of inductive coating diminish 
as G is reduced from 0.5 to 0.3. Then, of course, as G is reduced to zero, the corresponding 
curves would be identical to those shown in figures 2a and 2b. 

A capacitive-type coating corresponds to taking negative real values of G. Cross-section 
patterns are shown in figures 10c, lOd, and lOe for #=—0.1, —0.3, and —0.5, respectively. 
Here, it is interesting to note that the ripples are very small for both E- and //-plane patterns. 
At the same time, the forward scattering lobe is virtually unaffected by changes of G right 
through the range from +0.5 to —0.5. 

In figures 11a to lie, cross-section curves are shown to illustrate the effect of a conductive- 
type coating material such that Z=\Z\ exp {it/4) or 6r=(l— i)|6r|/-\/2 or <i>=45° (which are all 
equivalent statements). It is immediately evident that such a coating has the important 
effect of reducing the scattering cross sections for all directions except in the forward direction. 
Furthermore, the standing wave patterns are generally decreased by the presence of the coating. 
Comparisons between geometric optics and the exact theory for conductive-type coatings are 
shown in figures 12a and 12b for in the range from 90° to 180°. The agreement is reasonably 
good. 
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Figure 7a. 
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Figure 7b. 
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Figure 8a. 




Figure 8b. 



The effect of varying the phase angle <J> is illustrated for ka = 5 and |6?| = 0.10 in figures 
13a and 13b for jE'-plane scatter and in figures 14a and 14b for TJ-plane scatter. Those curves 
show quite clearly that the greatest reduction in scattering takes place when the surface im- 
pedance is purely real (i.e., $=0°), although marked reductions also occur when the coating 
is of a conductive type (i.e., $=45°). Similar-type curves are shown in figures 15a to 16b 
where &a=15, but otherwise the conditions are the same. 
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Figure 9a. 
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Figure 9b. 




Figure 10a. 




Figure 10b. 



6. Conclusions 

The numerical results presented graphically in this paper should provide some insight 
into the nature of scattering from coated spheres. Of principal importance is the fact that 
forward scattering is hardly influenced at all by the presence of reactive or lossy-type coatings. 
On the other hand, the magnitude of the backscattering cross sections is significantly modified 
by the magnitude and phase of the effective surface impedance of the coating. 
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Figure 10c. 




Figure lOd. 
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Figure lOe. 




Figure 11a. 



7. Appendix. Outline^of [Solution for Scattering of a Plane Wave From a 
Sphere With^Impedance Boundary Conditions 

The geometry of the problem is defined in terms of spherical coordinates (r, 9, </>) and 
common Cartesian coordinates (x, y, z) with common origin and common polar axis as indicated 
in figure 1. Thus, the surface of the sphere is at r=a and the incident plane wave is defined by 



E t x =E e- m =E e'- ikTCOBe y 

TTi __^2 p-Wz — ^O e -ikr coa 



(A-l) 
(A-2) 



for a time factor exp (iwt) and where 77 =120x. The total field can be written as a super- 
position of TM (transverse magnetic) and TE (transverse electric) waves. Following Harring- 
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ton [1961], the field components may be written in terms of potential functions A and F, 
respectively, for the two sets. Explicitly, 



a+ieoo\dr 2 ) 



H( ^ik(^ +k2 ) F ' 



H ( ( e) 



1_ b 2 A ; 

(a+itco)rbrbd 

1 bA 



$ r sin d<£ 

1 bF 
r sin 6 b<t> 

1 b 2 F 



Ef = 



ificor brbd 
1 



d 2 .4 



(<r+ito))r sin 9brb<j> 



H * ~ r bd 

** ~r bd 



Hi M =^ 



1 



b 2 F 



(A-3) 

(A-4) 
(A-5) 
(A-6) 
(A-7) 
(A-8) 
(A-9) 
(A-10) 
(A-ll) 
(A-12) 



(ilxo>)r sin Obrb<t> 
The total field in any region is the superposition of the TM and TE waves. For example, 

E t =El'»+E} n . (A-13) 




Figure lib. 




Figure 11c. 



310 



I 



1.0 



cr e (£) 



1.0 



i.o - 



0.1 



- 1 1 1 


1 1 1 1 1 r 


_ 


E-PLANE 




ka ■ 15 





GEOMETRIC OPTICS 




$ ■ 45° 


— 


IGI-- 0.15 


•/\/\. 


y ^^ / Ls^i-^~^- 


- 


— 


— 


i G 1= 0. 10 —zz 


- 


~ 


— 


|Ghb.05 


\-r- ~ ^7 




I ! ! 


1 1 1 1 1 



90° 100° 120° 140° 160° 180° 

9 (DEGREES) 



1.0 



1.0 — 



cr h (#) 



I 



— 1 1 


1 1 1 1 1 1 - 






_ 


H -PLANE 




ka = 15 


- 


GEOMETRIC OPTICS 




c£ -- 45° 





|G| = 0.15 


- 







|G| = 0.I0 


- 


|G|=0.05 




1 1 1 ! 1 1 



Figure 12a. 
The well-known addition theorem [Harrington, 1961], 
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Figure 12b. 



A 

i9 =Si"(2n+l)^P, (cose), 



permits us to write 



2?' r =cos0sin 02?i=2? o 



kr 

cos <t> d 
ikr dd 



( e -1*r oos •) ^ 



- iE ^! * ±i i-*(2n+l)J n (kr)P* (cos 6), 



(A-14) 

(A-15) 
(A-16) 



where the notation is that used by Schelkunoff [1943]. J n (kr) is a spherical Bessel function 
defined by 

J n (kr)=(^J 2 J n+m (kr), (A-17) 

where J n+U2 is a cylindrical Bessel function, Z J * is an associated Legendre polynomial denned 
here by 

P\ (cos B) = ~ P n (cos 6) , (A-18) 



in terms of the Legendre polynomial. The potential function A f , which gives rise to E r , is 
easily found to be 

(A-19) 



A*= -*k^± ± a„J n (kr)Pl (cos 6), 



M" B-l 
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Figure 13a. 
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Figure 13b. 



where 



dn = 



%- n (2n+\) 



n n(n+l) 



By a similar procedure, it is found that 



pi= _E^m^ ^ a J nikr)p i (cos e) (A _ 20) 

fC n-l 

For the region r>a, it is possible to write 

A=A f +A s , 
F=F i +F s , 

where the superscript s is associated with the scattered field. 

A 

A Clearly, A s and F s must contain the spherical Hankel functions H„ 2) (kr) in place of the 
J n (kr)'s, because the scattered field must behave as an outgoing wave at infinity. Thus, for 
the region r>a, 



and 



A=J2A n = - ^^* 2 [OnJ. (kr) + b n mf {kr) }P« (cos 6) , 



(A-21) 



F=J2 F^-^P S [flU.(*r)+«.fl? > (*r)]P. (0 (coa*), (A-22) 

7/ = 1 IC ft=l 



where b n and c n are coefficients unknown as yet. 
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Figure 14a. 
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Figure 14b. 



It is now convenient to introduce the concept of wave impedance such that relations 
between the tangential field components at r=a are to be prescribed. For example, for a 
TM mode, of order n. we insist that 



E^ = -Z^H^ and HJ^glg E? at r=a, 
while for a TE mode, of order n, 

H^=Y^E^ and fij»>= -^ flj« at r=a. 

These latter four conditions are automatically satisfied if 

_1 i>A n 

i 
and 



ie co dr 



--7J n) A n 



1 dF 
ipoca or 



(A-23) 



(A-24) 



(A-25) 



(A-26) 



Employing the boundary conditions on A n and F nf defined above, it is a simple matter to 
show that 



and 
where 

and 



_ i-»(2n+l) ff 

" 71(71+1) Dn * 



Cn = 



i(n+] 
i-»(2n+l) 

71(71+1) 



o _ J' n (ka)—iA n J n (ka) 



<?„=- 



fl»'(io)-»A.fl? , (4o) 

A A 

J' n (ka)—i8 H J n (ka) 



H™'(ka)-i8 n H™(ka) 

where fc=(e /to) l w, A B =Z (B 7»to, and 5„=i r( " ) 7;o. 
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Figure 15a. 
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Figure 15b. 



The solution of the first part of the problem has now been obtained. Namely, the fields 
may be expressed in terms of the normalized impedance A n and normalized admittance 8 n . 
Of particular interest are the far-zone scattered fields such that &r»l. It is then a simple 
matter to show that 



iE 



and 

where 

and 



El^+j^ e~ ikr cos 4>P(e)c±itoH% 



- (2ti+i) r Pi (cose) c d_ pi , e) l 



(A-31) 
(A-32) 
(A-33) 

(A-34) 



In order to arrive at an explicit solution, it is necessary to obtain formulas for the impedance 
Z {n) and the admittance Y (n) for TM and TE modes, respectively. This subject has been dis- 
cussed in an earlier paper by the author [Wait, 1963] and an equivalent problem has been treated 
by Garbacz [1962]. For purposes of the present paper, it is assumed that \Z {n) = 1/Y (n) =Z is\ 
specified complex constant independent of mode number n. By definition, Z becomes the sur- 
face impedance at the boundary r=a as discussed in the main text. 



The authors appreciate the extensive advice of K. P. Spies and L. A. Berry in connection 
with the computing work. In particular, the program for the spherical wave functions, as 
developed by L. A. Berry, was used extensively. 
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Figure 16a. 




Figure 16b. 



Note Regarding the Figures 

The vertical scale in a number of the figures is shifted in order to avoid overlapping of the curves when the 
ka value is changed. To prevent confusion, one should remember that <r(0) tends toward ira 2 when 0=180°. 
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